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ABSTRACT 

Using cosmological hydrodynamical simulations, we investigate the prospects of the thermal 
Sunyaev-Zel'dovich (tSZ) effect to detect the missing baryons in the local universe. We find that 
at least 80% of the tSZ luminosity is generated in collapsed structures, and that ~ 70% of the re- 
maining diffuse tSZ luminosity (i.e., ~ 15% of the total) comes from overdense regions with S gas >10, 
such as filaments and superclusters. The gas present in slightly overdense and underdense regions 
with S gas < 10, despite making up 50% of the total baryon budget, leaves very little tSZ signature: 
it gives rise to only ~ 5% of the total tSZ luminosity. Thus, future Cosmic Microwave Background 
(CMB) observations will be sensitive to, at best, one half of the missing baryons, improving the cur- 
rent observational status, but still leaving one half unobserved. Since most of the tSZ is generated in 
haloes, we find a tight correlation between gas pressure and galaxy number density. This allows us 
to predict the CMB Comptonization from existing galaxy surveys and to forecast the tSZ effect from 
the local structures probed by the Two Micron All Sky Survey (2MASS) galaxy catalog. 
Subject headings: cosmology: cosmic microwave background, observations - methods: numerical - 
galaxies: clusters: general 



1. INTRODUCTION 

Recent cos mological ob servations like galaxy surveys 
(e.g.,2 dFGRS, ICole et a.l I lpT05l) SPSS. iTegrnark et al.l 
(2004)), Supernovae (e.g., SNLS. lAstier et all (2006)) or 
Cosmic Microwave Backg round (CMB) (e.g., WMAP, 
iSpergel et al.1 l)2003[ 12006ft ) give strong support to the 
LCDM model of the universe, in which only ~ 18% of 
matter is in the form of baryons. Baryons are mostly 
probe d by CMB, Lym an-q forest (e.g., iCroft et alJ 
i|1998j) : iMcDonald et al.1 l|2005fl '). and X-rays observa- 
tions. The first two probes show that the Universe, at 
much younger epochs, contained an amount of baryons 
that exceeds by a factor of ~ 9 the baryo ns seen in 
today's universe, (jFukugita fc Peebles! 12004) . This is 
the missing ba r yons proble m. Indeed, according to 
ICen fc Ostrikerl |l999l 120051) . about half of the baryons 
should be in the form of a Warm Hot Intergalactic 
Medium (WHIM), at relatively high temperature (~ 
10 5 — 10 7 K) but low densities. Only a small frac- 
tion of the total gas is hot enough (>10 8 K) to en- 
able X-ray detection via bremmstrahlung: this mecha- 
nism led in the nineties to the identification of mas- 
sive galaxy clusters f r om R OSAT X-ray observations 
(e.g.. iBohringer ~ al. (2000J)), and more recently has 
motivated the search (and tentative detections) of dif- 
fuse X-ray emission from filaments and supercl usters , 
l|Kull fc Bohringerl Il999l iZannacosta et all 120021) . Al- 
though further pieces of evidence for local gas have been 
found in quasar spect ra, in the form of Oxyg en absorp- 
tion features at z ~ (jNicastro et al.l2005alH ). it is clear 
that we do not have a complete map of the baryon dis- 
tribution in the local Universe. 

With the advent of high resolution and high sensi- 
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tivity CMB experiments it has been proposed to trace 
the local gas using the therma l Suny aev-Zerdovich ef- 
fect ftSZ. lSimvarv fc ZeldovicH lff98fj) ). The tSZ effect 
describes the Compton scattering of CMB photons off 
hot electrons, in which a net transfer of energy from 
electrons to radiation distorts the Black Body spectrum 
of the CMB. The relative intensity fluctuations intro- 
duced in this scattering are redshift independent, and 
their amplitude is proportional to the integrated elec- 
tron pr essure along the line of sigh t. Observational ef- 
forts b y FCenova-Santos et alJ l|2005D and lBattistelli et alJ 
( 2006) have detected strong decrements in Corona Bore- 
alis, and interpreted them as Compton scattering being 
generated by extended gas in a fil ament aligned with 
the line of sight. Other studies lHansen et aT] 120051 
iDolag et all 120051: lAtrio-Barandela fc MiicketJ 12006(1 ad- 
dressed this problem from a theoretical point of view, 
trying to characterize the tSZ signal generated by non- 
collapsed baryons. 

Here we use a hydrodynamical numerical simulation to 
explore the prospects of future tSZ observations to search 
for the missing baryons. We find that, as the tSZ effect 
is sensitive to gas pressure, most of the tSZ luminosity 
(around 80%) is generated in collapsed structures, which 
host only ~30% of the total baryon mass. Most baryons 
are in environments with low overdensities, giving rise 
to very little tSZ signal: more than 60% of baryons are 
in regions with S gas < 20, which contribute only ~10% 
of the total tSZ luminosity. More than two thirds of 
the diffuse tSZ luminosity (not generated in collapsed 
structures) is coming from overdense environments with 
Sgas > 10 such as filaments or superclusters. The signifi- 
cant amount of baryons present in underdense regions (~ 
18% where 8 g a S < 0) is practically invisible to the tSZ ef- 
fect (~ 0.5 % contribution to tSZ luminosity). While the 
tSZ will be, at best, sensitive to <5 gas > 10 (i.e. superclus- 
ters and filaments), which corresponds to ~ 50% of all 
baryons, X-rays are only sensitive to the core centers of 
galaxy clusters in our neighborhood which corresponds to 
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Fig. 1. — (a) Relative contribution to the total tSZ luminosity of regions with gas density contrast bigger than <5 9as . The vertical solid 
line indicates the average gas overdensity in haloes, (b) Relative contribution to the total baryon mass of regions with gas density contrast 
bigger than <5 ga3 . The vertical solid line indicates the average gas overdensity in haloes. Note the asymmetry with respect to the previous 
(left hand side) panel, (c) Ratio of gas pressure produced in galaxy groups (small haloes) over pressure produced by diffuse gas and galaxy 
groups, versus total pressure (including galaxy clusters), averaged on scales of 12 h Mpc. The solid vertical line corresponds to regions 
with average matter density. In overdense regions, the pressure of haloes always dominates the diffuse pressure. 



less than 5% of all baryons, ijFukueita fc Peeblesl feOtM). 

Because in overdense regions most of the tSZ is gen- 
erated in haloes, we find a tight correlation between the 
gas pressure giving rise to the tSZ and the number galaxy 
density field, p e oc n 2 ^ 2 , . This allows us to build tSZ tem- 
plates and to make predictions for the tSZ signal from 
existing galaxy surveys. 

2. NUMERICAL SIMULATIONS 

We ran a cosmological simulation of a LCDM cosmol- 
ogy (Ct m = 0.3, A = 0.7, fl b = 0.045, h = 0.7, a 8 = 
0.9, n = 1) with a box of 200 Mpc/h side, 1024 3 hy- 
dro grid cells, and 512 3 dark m a tter p articles using the 
TVD+PM code of iTrac fc P en (20QJ). The comoving 
grid spacing is 195 kpc/h and the dark matter particle 
mass resolution is 5 x 10 9 M Q //i. The Eulerian hydro 
algorithm computes the spatial and temporal changes 
in the conserved mass, momentum, and total energy. 
The thermal energy is calculated by subtracting the ki- 
netic energy from the total energy, while the pressure 
and temperature are related to the thermal energy by 
the standard equation of state. The simulation is adia- 
batic. We recall the definition of the the Comptonization 
parameter (y = J drax n e (kBT e )/(m e c 2 ), with o~t the 
Thomson cross-section, ks the Boltzmann constant and 
m e ,n e ,T e the electron mass, density and temperature, 
respectively), in order to introduce the derivative dy/dr 
as a proxy for electron pressure, (dy/dr cx p e oc n e T e ). 
Likewise, we define Y = J d 3 r dy/dr, as a proxy for tSZ 
luminosity, since both quantities are also proportional. 
Bound dark matter particles classified as halos are iden- 
tified using an ellipsoidal overdensity threshold of 200 
times the average density. This defines the halo bound- 
aries which are used to distinguish between tSZ gener- 
ated in collapsed structures and in the diffuse phase. 

3. WHERE IS THE TSZ PRODUCED? 

We distinguish tSZ luminosity being generated in cells 
belonging to clusters of galaxies (M c i > 5xl0 13 h~ 1 M Q ), 
to small haloes (also referred to as galaxy groups, and de- 



fined as all haloes resolved in the simulation with masses 
below 5 x 10 13 h~ 1 M ( ? ) ), and to a diffuse gaseous phase 
(defined as all gas cells not belonging to any halo) . Given 
the dark mass particle resolution in our simulation, we 
choose our mass threshold for the halo definition to be 
~ 10 12 h~ 1 M Q : all haloes below this limit are regarded 
as diffuse gas, and hence our estimates for the diffuse tSZ 
contribution should be regarded as optimistic. 

Figure (JTJi) shows the cumulative distribution of the 
pressure in the box versus gas density contrast. Col- 
lapsed structures show a halo-mass- weighted average gas 
overdensity of ~ 93, (marked by the vertical solid line). 
When integrating -i.e., volume weighting- the tSZ lumi- 
nosity in the three cell subsets defined above (galaxy clus- 
ters, small haloes and diffuse gas), we find that ~ 70% 
of total tSZ luminosity is generated in galaxy clusters, 
and, out of the ~ 30% remaining, around one third of 
it (i.e., ~ 10% of the total) is generated in small haloes. 
This leaves the diffuse phase a ~ 20% contribution, even 
though it hosts ^70% of the total baryonic mass. Fur- 
ther, one fourth of this diffuse gas is located in under- 
dense regions, whose tSZ luminosity is negligible (< 1% 
) . This asymmetry of the pressure and baryon mass dis- 
tributions versus gas overdensity is explicitely shown in 
Figures |^i,b). 

This suggests that the detection, via the tSZ effect, 
of the diffuse gas will be hampered by the presence of 
small (and plausibly unresolved) haloes: this is shown in 
Figure Jflj), where the ratio of the gas pressure gener- 
ated in small haloes over the sum of the contributions 
from the diffuse gas and small haloes is plotted versus 
the total gas pressure, (in all cases the pressure is com- 
puted in cells of 12 /i _1 Mpc side). The vertical line marks 
the total pressure corresponding to the average galaxy 
number density. We see that, in slightly overdense re- 
gions hosting large pressure (right of the vertical line), 
the tSZ luminosity outside galaxy clusters is preferen- 
tially generated in smaller haloes rather than in diffuse 
gas, i.e., when searching for the tSZ signature of diffuse 
gas in overdense regions one must carefully subtract the 
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Fig. 2. — (a) Pressure — galaxy number density correlation, according to SSHJ (filled circles and solid line) and K04 (filled triangles, 
dashed line) prescription to populate haloes with galaxies, (b) Correlation between Comptonization parameter y and projected galaxy 
overdensity (K04 prescription) obtained after projecting a simulated sphere of of 300 h~ 1 Mpc radius, (c) Rayleigh- Jeans (RJ) expectation 
for our prediction of the local tSZ power spectrum computed from 2MASS galaxy catalog and the y-galaxy correlation (filled circles). The 
tSZ power spectrum (RJ) from our simulated volume is given by the solid line, and decomposed into the halo (dashed line) and diffuse 
(dot-dashed line) contributions. 



contribution from small haloes. Projection effects in fu- 
ture CMB maps will further enhance this contamination. 
Only when considering the very overdense regions of the 
universe (far right end of this plot), the relative weight of 
the diffuse phase takes over, since in such environments 
there are no small haloes (all haloes are clusters). Note 
however that the mass threshold for our cluster definition 
(5 x 10 13 }i~ 1 Mq) corresponds to abundant and not-so- 
tSZ-bright objects, and that these sources may still be 
responsible for a significant amount of confusion noise. 

4. THE PRESSURE - HALO CORRELATION 

Since more than 80% of the tSZ luminosity is generated 
in collapsed structures, there must be a correlation be- 
tween the gas pressure and the number density of haloes, 
or the gas pressure and the galaxy number density. 

To investigate this correlation, we populate the haloes 
in th e simulation using two dif ferent algorithms: the first 
one ijScoccimarro et ahl ()2001[) . hereafter SSHJ) charac- 
terizes the halo galaxy populat ion by a binom i al dis - 
tribution, whereas the second l)Kravtsov et alJ ^2004^1 . 
hereafter K04) us es a Poissonian distribution. As in 
Smit h~et alJ (|2P06), both algorithms are normalized to 
give the same average galaxy number density, {(n g ) ~ 
3 x 1CP 2 /i 3 Mpc~ 3 ), and this requires slightly different 
choices for the halo minimum mass hosting a galaxy. 
As shown in Figure ((30, the correlation found between 
dy/dr and galaxy number density (computed in cells of 
12/i -1 Mpc side) is very similar in both cases, and well 
fitted by a power law of index n ss 2.1 — 2.3 for the SSHJ 
and K04 methods (filled circles and filled triangles, solid 
and dashed lines, respectively). The slope obtained in 
both cases is close to the prediction that a polytropic 
gas model {p cx p 1 ) provides for a self-gravitating sys- 
tem, (7 = 2). 

As current all sky maps of the local galaxy distribu- 
tion do not have distance information, we want to inves- 
tigate the correlation between projected galaxy density 
and y. We extend our simulated volume to a box of 
600 h _1 Mpc side, exploiting the fact that the simulation 



box has periodic boundary conditions. We then place 
an observer in the center and project on a HEALPix 
lIGorski et alJ f2005) map of resolution parameter N^e 
= 512 a sphere of 300 h _1 Mpc radius. We produce maps 
of y and galaxies, smooth them with a Gaussian beam of 
FWHM = 12.6 arcmins, and sort the pixels in increas- 
ing projected number of galaxies. After binning them 
in groups of 32 pixels (for display purposes), we obtain 
a correlation between y and the projected galaxy over- 
density (n g (n)/{rin), n denoting a direction on the sky), 
shown in Figure pfc). 

We find that for slight galaxy overdensities (few times 
the background density), there is a tight correlation be- 
tween galaxy number density and the Comptonization 
parameter y. This suggests using existing galaxy cata- 
logs to predict the degree of Comptonization of the CMB 
sky. The volume used in our projection is close to that 
sampled by the Two Micron All Sky Survey (hereafter 
2MASS, Uarrett et all l)2003j) L We thus use a map of 
projected galaxies from this s urvey (constructed as in 
Hcrnandc z-Monteagudo et alJ {2004)) to make a predic- 
tion of the y sky. When comparing the galaxy maps 
from 2MASS and our simulation we find that they both 
show similar galaxy densities and clustering properties. 
We approximate the correlation given in Figure (03 ) by 
five power laws (given by the solid lines), and invert the 
2MASS-based galaxy catalog into a y map. 

The overall y normalization is set by imposing that 
the 256 highest density pixels have a tSZ decre- 
ment of -73 ± 17 pK, as found at 94 GHz by 
iHernandez-Monteagudo et alJ i|2004D (all maps have 
been convolved by the same Gaussian beam) . The power 
spectrum of the resulting y map (in units of [/iK] 2 at 
Rayleigh- Jeans frequencies) is shown in Figure (J2t) by 
filled circles: its amplitude and shape is very close, at 
high multipoles, to the power spectrum of the y map 
produced from our simulation (solid line). They only 
differ clearly at large scales (£ < 10), for which the y 
power spectrum estimated from 2MASS is higher. This 
prediction for the local tSZ power spectrum is remark- 
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Fig. 3. — (a) All-sky diffuse tSZ contribution in our simulated volume, as oppossed to the total tSZ contribution, (panel (b)). 



ab ly close both in shape and amplitud e to th ose obtained 
bv lHansen et all 1)2005(1 : iDolae et all l|2005[) via a inde- 
pendent approach. Our computation should be a valid 
prediction for the tSZ power at low is, but at high is the 
neglected high-rcdshift cluster population will dominate. 

5. PROSPECTS FOR DIFFUSE TSZ DETECTION 

Since we have separated in our simulation the gas be- 
longing to collapsed haloes from the gas in a diffuse 
phase, it is possible to compute y maps and power spec- 
tra of the halo and the diffuse components separately. 
The dot-dashed line at the bottom of Figure J2t) pro- 
vides the power spectrum of the diffuse gas in our sim- 
ulation, whereas the dashed line corresponds to the tSZ 
power spectrum generated by gas located in haloes (note 
its proximity to the total contribution [solid line]). A 
map of diffuse tSZ is given in Figure (JHtO, and com- 
paring to the total tSZ contribution (Figure (03)). We 
conclude that most of the diffuse tSZ signal visible on 
the CMB sky (and more than 90% of the total tSZ lu- 
minosity) is generated in overdense regions traced by 
the clusters and superclusters, (i.e., the diffuse warm 
gaseous component present in voids leaves negligible tSZ 
signature). The amplitude of the diffuse phase is at 
best ten times (a hundred times in CYs) smaller than 
the tSZ signal generated in haloes, even at large angu- 
lar scales. In particular since galaxies trace superclus- 
ters and filaments (corresponding to S gas > 10), the 
diffuse gas in these regions could be targeted by tSZ 
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observations with specifications similar to those of e.g. 
ACT (http://www.hep.upenn.edu/act/ ). These re- 
gions contain 50% of the baryons: thus tSZ observations 
may decrease the ratio of unseen baryons in the local 
universe from ~ 9 to ~ 2. In practice, the limiting factor 
in detecting this signal will be confusion noise from un- 
resolved haloes. Note that the ratio of tSZ luminosities 
of the collapsed and diffuse baryon components does not 
equal the ratio of the corresponding tSZ-induccd angular 
power spectra, since the latter quantities depend on the 
contrast of tSZ sources against the CMB background. 
This makes the search for the (diffuse) missing baryons 
even more challenging. The detection of part of the dif- 
fuse tSZ may be feasible in nearby superclusters after 
masking out all compact sources and well characterized 
nearby galaxy clusters. However, we must conclude that 
the tSZ, despite of providing a new tool to unveil the 
presence of unseen warm and hot gas, will still miss a 
significant fraction of the baryons in the local Universe. 



We thank R.E.Smith and R.Sheth for useful dis- 
cussions. We a cknowledge the use of HEALPix 
l|G6rski et alJ 12005]) package, and the Two Micron All 
Sky Survey. CHM is supported by NASA grants ADP03- 
0092 and ADP04-0093 and NSF grant PIRE-0507768. 
RJ is partially supported by NSF grant PIRE-0507768 
and by NASA grant NNG05GG01G. LV is supported by 
NASA grants ADP03-0092 and ADP04-0093. 



Hansen, F. K., Branchini, E., Mazzotta, P., Cabella, P., & Dolag, 
K. 2005, MNRAS, 361, 753 

Hernandez-Monteagudo, C, Genova-Santos, R., & Atrio- 
Barandela, F. 2004, ApJL, 613, L89 

Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., & Huchra, 
J. P. 2003, AJ, 125, 525 

Kravtsov, A. V., Berlind, A. A., Wechsler, R. H., Klypin, A. A., 
Gottlober, S., Allgood, B., & Primack, J. R. 2004, ApJ, 609, 35 

Kull, A., Bohringer, H. 1999, A&A, 341, 23 

McDonald, P., et al. 2005, ApJ, 635, 761 

Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005. 

ApJ, 623, L63 
Nicastro, F., et al. 2005, ApJ, 629, 700 
Nicastro, F., et al. 2005, Nature, 433, 495 

Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001, ApJ, 546, 
20 

Sheth, R. K., & Tormen, G. 1999, Mon. Not. Roy. Astron. Soc, 
308, 119 



5 



Smith, R. E., Watts, P. I. R., & Sheth, R. K. 2006, MNRAS, 365, 
214 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

S nergel. D. N.. et al. 2006, ArXiv Astrophysics c-prints, 

arXiv:astro-ph/0603449 
Sunyaev, R. A., & z,eiaovich, Y. B. 1980, ARA&A, 18, 537 
Tegmark, M., et al. 2004, ApJ, 606, 702 



Trac, H., & Pen, U.-L. 2004, New Astronomy, 9, 443 
Zappacosta, L., Mannucci, P., Maiolino, R., Gilli, R., Ferrara, A., 

Finoguenov, A., Nagar, N. M., & Axon, D. J. 2002, A&A, 394, 

7 



